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Abstract 

I review the recent calculations and current status of the hadronic 
light-by-light scattering contribution to muon g-2. In particular, I 
discuss the main results obtained in a recent work together with Ed- 
uardo de Rafael and Arkady Vainshtein where we came to the estimate 
a^ LbL = (10.5=L2.6) x 10~ 10 . How the two-photon physics program of 
low energy facilities can help to reduce the present model dependence 
is also emphasized. 



1 Invited talk at "Photon09, International Conference on the Structure and the Inter- 
actions of the Photon", May 11-15 2009, DESY, Hamburg, Germany. 



1 Introduction 



One momenta configuration out of the six possible ones contributing to the 
hadronic light-by-light to muon g-2 is depicted in Fig. [1] and described by 
the vertex function 
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where q — > is the momentum of the photon that couples to the external 
magnetic source, q = p 2 — Pi = —k\ — k 2 — k% and m is the muon mass. 




Figure 1: Hadronic light-by-light scattering contribution. 

The dominant contribution to the hadronic four-point function 

W^(qMMM) = 
t 3 [ d 4 x [ d 4 y [ d 4 z e K-ki-x+k 3 . y +h 2 .z) ^ | T [ V »(0)V u (x)V p (y)V u (z)} |0> (2) 



comes from the three light quark (q = u, d, s) components in the electromag- 
netic current = qQj^q (x) where Q = diag(2, — 1, — 1)/3 denotes 
the quark electric charge matrix. We are interested in the limit q — > where 
current conservation implies 
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Therefore, the muon anomaly can then be extracted as 
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Here I discuss the results of [1] and [2]. Previous work can be found in 
PUMIHEiEidniimil^and recent reviews are in [El El ESI [16] . 

The hadronic four-point function n^ p<T (g, hi, k 3 , k 2 ) is an extremely dif- 
ficult object involving many scales and no full first principle calculation of it 
has been reported yet -even in the simpler large numbers of colors N c limit of 
QCD. Notice that we need it with momenta fci, k 2 and k 3 varying from to 
oo. Unfortunately, unlike the hadronic vacuum polarization, there is neither 
a direct connection of a HLbL to a measurable quantity. Two lattice groups 
have started exploratory calculations [TTJ, [18] but the final uncertainty that 
these calculations can reach is not clear yet. 

Attending to a combined large number of colors of QCD and chiral per- 
turbation theory (CHPT) counting, one can distinguish four types of con- 
tributions [19]. Notice that the CHPT counting is only for organization of 
the contributions and refers to the lowest order term contributing in each 
case. In fact, Ref. pQ shows that there are chiral enhancement factors that 
demand more than Nambu-Goldstone bosons in the CHPT expansion in the 
light-by-light contribution to the muon anomaly. See more comments on this 
afterwards. 

The four different types of contributions are: 

• Nambu-Goldstone boson exchanges contribution are 0(N C ) and start 
at 0(p 6 ) in CHPT. 

• One-meson irreducible vertex contribution and non-Goldstone boson 
exchanges contribute also at O(N ) but start contributing at 0(p 8 ) in 
CHPT. 

• One-loop of Goldstone bosons contribution are 0(1/N C ) and start at 
C(p 4 ) in CHPT. 

• One-loop of non-Goldstone boson contributions which are 0(1/N C ) but 
start contributing at 0(p 8 ) in CHPT. 
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Based on the counting above there are two full calculations [21 IH E] and 
[5J [7]. There is also a detailed study of the n° exchange contribution [S] 
putting emphasis in obtaining analytical expressions for this part. 

Recently, two new calculations of the pion exchange using also the orga- 
nization above have been made. In Ref. [10], the pion pole term exchange is 
evaluated within an effective chiral model. These authors also study the box 
diagram one-meson irreducible vertex contribution. The results are numeri- 
cally very similar to the ones found in the literature as can be seen in Table 
[TJ In Ref. [11] , the author uses a large N c model 7r°7*7* form factor with the 
pion also of-shell. This has to be considered as a first step and more work 
has to be done in order to have the full light-by-light within this approach. 
In particular, it would be very interesting to calculate the contribution of 
one-meson irreducible vertex contribution within this model. 

Using operator product expansion (OPE) in QCD, the authors of [12] 
pointed out a new short-distance constraint of the reduced full four-point 
Green function 

(0|T [V^k^V^V^-ih + k 2 + q))] | 7 ( 9 )> (5) 

when q — > and in the special momenta configuration — k\ ~ —k\ » 
—(ki + ks) 2 Euclidean and large. In that kinematical region, 

with k = (k\ — k 3 )/2 ~ k\ ~ — k ?J . See also [20J. This short distance 
constraint was not explicitly imposed in previous to [12] calculations. 



2 Leading in 1/N C Results 

Using effective field theory techniques, the authors of [9] shown that leading 
large iV c contribution to a HLbL contains an enhanced term at low energy by 
log 2 (M p /m n ) where the rho mass M p acts as an ultraviolet scale and the pion 
mass m n provides the infrared scale. 

„=*V) = f-) 3 k - + O (in ^) + 0(1)1 (7) 

\7r/ 487i - Z J* I m n V m n J 

This leading logarithm is generated by the Goldstone boson exchange con- 
tributions and is fixed by the Wess-Zumino-Witten (WZW) vertex 7r°77. In 
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Figure 2: A generic meson exchange contribution to the hadronic light-by- 
light part of the muon g — 2. 

the chiral limit where quark masses are neglected and at large N c , the coeffi- 
cient of this double logarithm is model independent and has been calculated 
and shown to be positive in [S]. All the calculations we discuss here agree 
with these leading behaviour and its coefficient including the sign. A global 
sign mistake in the 7r° exchange in [31 HJ Ej was found by [H [9] and confirmed 
by [6l[7] and by others [HI [22]. The subleading ultraviolet scale /z-dependent 
terms [9], namely, \og{^/m n ) and a non-logarithmic term are model 

dependent and calculations of them are implicit in the results presented in 
[3, HI EJ H2]. In particular, k([jl) contains the large N c contributions from 
one-meson irreducible vertex and non-Goldstone boson exchanges. In the 
next section we review the recent model calculations of the full leading in 
the 1/N C expansion contributions. 

2.1 Model Calculations 

The pseudo-scalar exchange is the dominant numerical contribution and was 
saturated in [3], SI El El [TJ, [HI [TQl [11] by Nambu-Goldstone boson's exchange. 
This contribution is depicted in Fig. 2 with M = 7r°, 77, 77'. The relevant 
four-point function was obtained in terms of the off-shell 7r°7*(fci)7*(fc3) form 
factor ^(kf, fcf) and the off-shell vr°7*(fc2)7(g = 0) form factor 0) mod- 

ulating each one of the two WZW 7r°77 vertex. 

In all cases several short- distance QCD constraints were imposed on these 
form-factors. In particular, they all have the correct QCD short- distance 
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Table 1: Results for the 7r°, 77 and 7/ exchange contributions. 



Reference 


10 


10 x a 




7T° only 


7T°, 77 and 7/ 


ME] 


5.7 


8.3 ± 0.6 


mm 


5.6 


8.5 ± 1.3 


[8] with h 2 = 


5.8 


8.3 ± 1.2 


[8] with /i 2 = -10 GeV 2 


6.3 




m 


6.3 ~ 6.7 




m 


7.2 


9.9 ± 1.6 


m\ 


7.65 


11.4±1.0 



behaviour 

HQ 2 ,Q 2 )-^~ 2 and ^(Q 2 ,0)^^ (8) 

when Q 2 is Euclidean and large and are in agreement with 7r°7*7 low-energy 
data IE They differ slightly in shape due to the different model assumptions 
(VMD, ENJL, Large N c , NxQM) but they produce small numerical differ- 
ences always compatible within quoted uncertainty ~ 1.3 x 10~ 10 -see Table 

m 

Within the models used in P H El El El QUI [EE] , to get the full contri- 
bution at leading in 1/N C one needs to add the one- meson irreducible vertex 
contribution and the non-Goldstone boson exchanges. In particular, below 
some scale A, the one-meson irreducible vertex contribution was identified in 
[HI [7J with the ENJL quark box contribution with four dressed photon legs. 
While to mimic the contribution of short-distance QCD quarks above A, a 
loop of bare massive heavy quark with mass A and QCD vertices was used. 
The results are in Table [2] where one can see a very nice stability region when 
A is in the interval [0.7, 4.0] GeV. Similar results for the quark loop below A 
were obtained in [31 H] though these authors didn't discuss the short- distance 
long-distance matching. 

In [51 [7] , non-Goldstone boson exchanges were saturated by the hadrons 
appearing in the model, i.e. the lowest scalar and pseudo-vector hadrons. 
Both states in nonet-symmetry -this symmetry is exact in the large N c limit. 
Within the ENJL model, the one-meson irreducible vertex contribution is 

2 See however the new measurement of the 77* — > ttq transition form factor by BaBar 
[23] at energies between 4 and 40 GeV 2 
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Table 2: Sum of the short- and long-distance quark loop contributions [5] as 
a function of the matching scale A. 



A [GeV] 


0.7 


1.0 


2.0 


4.0 


10 10 x a HLbL 


2.2 


2.0 


1.9 


2.0 



Table 3: Results for the axial-vector exchange contributions from [31 HI [6] 
and [5HU. 



References 


10 10 x a HLbL 


ME] 

mm 


0.17 ± 0.10 
0.25 ± 0.10 



related trough Ward identities to the scalar exchange which we discuss below 
and both have to be included [5J [7J . The result of the scalar exchange obtained 
in m is 

a HLbL (Scalar) = -(0.7 ± 0.2) x lO" 10 . (9) 

The scalar exchange was not included in 0, HI El [8] . The result of the axial- 
vector exchanges in HI E] and m m can be found in Table El 

Melnikov and Vainshtein used a model that saturates the hadronic four- 
point function in (|2j) at leading order in the 1/N C expansion by the exchange 
of the Nambu-Goldstone ir°, rj, rj and the lowest axial-vector fx states. In that 
model, the new OPE constraint of the reduced four-point function found in 
[T2] mentioned above, forces the 7r°7*(g)7(j>3 = 0) vertex to be point-like 
rather than including a T{q 2 , 0) form factor. There are also OPE constraints 
for other momenta regions [21] which are not satisfied by the model in [T2] 
though they argued that this made only a small numerical difference of the 
order of 0.05 x 10 -10 . In fact, within the large N c framework, it has been 
shown [25] that in general for other than two-point functions, to satisfy fully 
the QCD short-distance properties requires the inclusion of an infinite num- 
ber of narrow states. 

3 Next-to- leading in 1/N C Results 

For the next-to- leading in 1/N C contributions to the a HLbL there is no model 
independent result at present and is possibly the most difficult component. 
Charged pion and kaon loops saturated this contribution in [31 HI El El CTJ • To 
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— Like 



Figure 3: Goldstone boson exchange in the model in [12] contributing to the 
hadronic light-by-light. 

dress the photon interacting with pions, a particular Hidden Gauge Symme- 
try (HGS) model was used in [31 HI E] while a full VMD was used in [HJ [7J. 
The results obtained are -(0.45±0.85) x 1(T 10 in [3] and -(1.9±0.5) x 10" 10 
in [5] while using a point-like vertex one gets —4.6 x 10~ 10 . Both models sat- 
isfy the known constraints though start differing at 0(p 6 ) in CHPT. Some 
studies of the cut-off dependence of the pion loop using the full VMD model 
was done in [5] and showed that their final number comes from fairly low 
energies where the model dependence should be smaller. The authors of [T2] 
analyzed the model used in [31 Hj and showed that there is a large cancel- 
lation between the first three terms of an expansion in powers of (m^/Mp) 2 
and with large higher order corrections when expanded in CHPT orders but 
the same applies to the 7r° exchange as can be seen from Table 6 in the first 
reference in [2j by comparing the WZW column with the others. The au- 
thors of [12] took (0 ± 1) x 10 -10 as a guess estimate of the total NLO in 
1/N C contribution. This seems too simply and certainly with underestimated 
uncertainty. 

4 Comparing Different Calculations 

The comparison of individual contributions in [31 HI El El El HOI ECU U2] has 
to be done with care because they come from different model assumptions 
to construct the full relevant four-point function. In fact, the authors of [10] 
have shown that their constituent quark loop provides the correct asymp- 
totics and in particular the new OPE found in [12J. It has more sense to 
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Table 4: Results for the full hadronic light-by-light contribution to a 



Full Hadronic Light-by-Light 


10 10 x a M 


[31 HIE] 


8.9± 1.7 


0121 


8.9 ± 3.2 


m 


13.6 ± 2.5 



compare results for a either at leading order or at next-to-leading order 
in the 1/N C expansion. 

The results for the final hadronic light-by- light contribution to a HLbL 
quoted in [HJ HI El El El ^] are i n Table HI The apparent agreement be- 
tween [31 HI E] and [51 E] hides non-negligible differences which numerically 
almost compensate between the quark-loop and charged pion and [12] are 
in Table HI Notice also that [31 HI E] didn't include the scalar exchange. 
Comparing the results of [HI E] and [12] , as discussed above, we have found 
several differences of order 1.5 x 10~ 10 which are not related to the new 
short- distance constraint used in [12]. The different axial- vector mass mix- 
ing accounts for —1.5 x 10~ 10 , the absence of the scalar exchange in [T2] 
accounts for —0.7 x 10~ 10 and the use of a vanishing NLO in 1/N C contri- 
bution in [12] accounts for —1.9 x 10~ 10 . These model dependent differences 
add up to —4.1 x 10~ 10 out of the final —5.3 x 10~ 10 difference between [HI E] 
and [12] . Clearly, the new OPE constraint used in [12] accounts only for a 
small part of the large numerical final difference. 



5 Conclusions and Prospects 

To give a result at present for the hadronic light-by-light contribution to 
the muon anomalous magnetic moment, the authors of [1], from the above 
considerations, concluded that it is fair to proceed as follows 
Contribution to a HLbL from tt° , rj and rf exchanges 

Because of the effect of the OPE constraint discussed above, we suggested to 
take as central value the result of Ref. [12] with, however, the largest error 
quoted in Refs. [HI E] : 

a HLbL (vr ,ij,rf) = (11.4 ± 1.3) x 1(T 10 . (10) 

Recall that this central value is quite close to the one in the ENJL model 
when the short-distance quark loop contribution is added there. 
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Contribution to a HLbL from pseudo-vector exchanges 

The analysis made in Ref. [12] suggests that the errors in the first and second 
entries of Table 2 are likely to be underestimates. Raising their ±0.10 errors 
to ±1 puts the three numbers in agreement within one sigma. We suggested 
then as the best estimate at present 

a HLbL (pseudo - vectors) = (1.5 ± 1) x 10~ 10 . (11) 

Contribution to a HLbL from scalar exchanges 

The ENJL-model should give a good estimate for these contributions. We 
kept, therefore, the result of Ref. [5J [7] with, however, a larger error which 
covers the effect of other unaccounted meson exchanges, 

a HLbL (scalars) = -(0.7 ± 0.7) x 10~ 10 . (12) 

Contribution to a HLbL from a dressed pion loop 

Because of the instability of the results for the charged pion loop and unac- 
counted loops of other mesons, we suggested using the central value of the 
ENJL result but wit a larger error: 

a HLbL (7r-dressed loop) = -(1.9 ± 1.9) x 10~ 10 . (13) 

From these considerations, adding the errors in quadrature, as well as the 
small charm contribution 0.23 x 10~ 10 , we get 

a HLbL = (10.5 ±2.6) x 10~ 10 , (14) 

as our final estimate. 

The proposed new g^ — 2 experiment accuracy goal of 1.6 x 10~ 10 calls for 
a considerable improvement in the present calculations. The use of further 
theoretical and experimental constraints could result in reaching such accu- 
racy soon enough. In particular, imposing as many as possible short-distance 
QCD constraints [31 lU El El El [11] has result in a better understanding of 
the numerically dominant 7r° exchange. At present, none of the light-by-light 
hadronic parametrization satisfy fully all short distance QCD constraints. In 
particular, this requires the inclusion of infinite number of narrow states for 
other than two-point functions and two-point functions with soft insertions 
[25] . A numerical dominance of certain momenta configuration can help to 
minimize the effects of short distance QCD constraints not satisfied, as in 
the model in [T2] . 
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More experimental information on the decays it — > 77*, 7r° — > 7*7* and 
7t° — > e + e~ (with radiative corrections included EZ] ) can also help to 

confirm some of the neutral pion exchange results. A better understanding 
of other smaller contributions but with comparable uncertainties needs both 
more theoretical work and experimental information. This refers in particular 
to pseudo- vector exchanges. Experimental data on radiative decays and two- 
photon production of these and other C-even resonances can be useful in that 
respect. 

New approaches to the pion dressed loop contribution, together with ex- 
perimental information on the vertex 7r + 7r~7*7* in the intermediate energy 
region (0.5 — 1.5 GeV) would also be very welcome. Measurements of two- 
photon processes like e + e~ — > e + e~7r + 7r~ can be useful to give information on 
that vertex and again could reduce the model dependence. The two-gamma 
physics program low energy facilities like the experiment KLOE-2 at DA$NE 
will be very useful and well suited in the processes mentioned above which 
information can help to decrease the present model dependence of a^ LbL . 
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